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Abstract 

The flux of turbulent kinetic energy from large to small spatial scales is measured in a small 
domain B of varying size R. The pdf of the flux is obtained using a time-local version of Kolmogorov 
four-fifths law. The measurements, made at moderate Reynolds number, show frequent events 
where the flux is backscattered from small to large scales, their frequency increasing as R is 
decreased. The observations are corroborated by a numerical simulation based on the motion 
of many particles and on an explicit form of the eddy damping. 

PACS numbers: 47.27.Eq, 47.27.Gs 
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I. INTRODUCTION 



There is considerable interest in the global (spatially averaged) fluctuations in systems 
driven far from equilibrium QQSQ.of which ffmd turbulence prov.des a striking example 
[5|,|6|]. An essential aspect of three-dimensional turbulence is the cascade of energy from large 
to small scales, followed by dissipation at the smallest scales. Characterizing the energy flux 
is particularly important for turbulence modeling. The local rate of energy dissipation is 
known to fluctuate wildly pj. This work investigates flux of energy from large to small 
scales, averaged over a local region of finite extent. 

The spatially averaged value of the energy flux e over a subvolume of the fluid B of 
typical size R, is simply given by the rate of dissipated kinetic energy if the system is in a 
steady state. In that case the flux is necessarily positive and is directed from large to small 
scales. However, temporal fluctuations in this flux can be very substantial. In fact, it is has 
been demonstrated several times that energy may backscatter from small to large scales, 
leading to a negative value of the energy flux [ll . It is natural to expect that this effect 
should depend on the scale of the subsystem investigated. One of the aims of this work is 
to quantify the fluctuations of the energy flux measured over a subdomain of the flow, and 
in particular, its dependence on the subdomain's size. 

Under conditions of local isotropy, the ensemble averaged energy dissipation rate, e, is 
related to the third moment of the longitudinal velocity differences at a given scale, r : 

<(Au L (r)) 3 > = ~er (1) 



The recent theoretical 19 
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111 ] and numerical work [12j leads to a generalisation of Eq. 
which is local in space and time. More precisely, if one considers any finite subdomain of the 
flow, B, the average of (Aii/^r)) 3 over all directions of the vector r, and over all points in 
B, is equal to z g-£ j g(t)r, where Eb is the energy flux towards small scales in the subdomain 
B : 

1 d(v*) ((Au L (r,t)f) B ,n 

£B{t) = 2^r - — => — (2) 

Here (•••)b,q denotes the average over the subvolume B, and over all possible directions 
(Q). The derivation of Eq. El is based on a rigorous energy balance for weak solutions of 
the Navier-Stokes equations, in the Re — > oo limit jl(3]. In this sense the quantity is 



really the instantaneous, inertial range dissipation rate, interpreted here as an energy flux 
(defined as positive towards small scales). Eq. El justifies the intuitive estimates of the rate 
of change of energy in a volume, by a straightforward averaging with the help of Eq. Q even 
when the key assumption of local isotropy is not satisfied. The spherical avera ge i n Eq. |2] 
permits the recovery of the isotropic sector from an arbitrary (anisotropic) flow |l2|. 

In practice, the average over all directions is replaced here by an average over particle 
pairs along all azimuthal directions in an instantaneous snapshot within the subdomain B 
of interest. The magnitude and sign of the numerator on the RHS fluctuates in time, though 
it must be negative on the average if the system is in the steady state. This assures that 
the time average of the flux is positive. 

Averaging over all pairs of particles in a sub-domain of the fluid, as done here experimen- 
tally, has also been proposed as a way of estimating the energy dissipation in a numerical 
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. In such an ap- 



scheme, based on the "smooth particle hydrodynamics" (SPH) method 
proach the values of various hydrodynamic fields at a point x are obtained by interpolating, 
with the help of all the neighbouring particles, within a smoothing distance h from the point 
x. An expression for the eddy-damping, very similar to Eq. El has been postulated in ^| 
(see also discussion below). Implementing numerically this eddy-damping in a SPH code, 
as done in ^| leads to a Large Eddy Simulation scheme, based on particles moving with 
the flow (lagrangian particles), where the smoothing length h is the smallest scale resolved 
numerically. The present experimental setup allows an explicit measurement of the en ergy 
flux from Eq. El and a comparison with the postulated form of the eddy-damping |l4| . 
The experimental results demonstrate that the postulated form of the energy-dissipation 
term is very strongly correlated with the energy flux obtained from Eq. Furthermore, the 
numerical scheme based on SPH allows one to study the energy flux in a subdomain of the 
system as a function of size. The qualitative agreement between the experimental and the 
numerical studies is excellent. 

II. EXPERIMENTAL METHOD 



The experimental arrangement is identical to the one used previously for studying surface 
turbulence and is shown in Fig. ^ Measurements are made 4 cm below the surface 

in a square tank of side length 1 m, filled with water to a height of 30 cm. The fluid is seeded 
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TABLE I: Turbulent quantities of interest measured in the experiment. 



Parameter 


Expression 


Measured Value 


Taylor Microscale A (cm) 


/ u rms 

\ 1 1 rlii /r)x\ 2 \ 
y \(OU x /UX) j 


0.32 


Taylor Microscale Reynolds Number Re\ 


V 


79 


Integral Scale Iq (cm) 


[ A x l {(u L (x+r)u L (x))) x 


3.5 


Large Eddy Turnover Time t; (s) 


lo 

Urms 


1.4 


Dissipation Rate Sdiss (cm 2 /s 3 ) 


15^(tf) 2 > 


8.3 


Kolmogorov Scale rj (cm) 


/ u a \l/4 


0.02 


RMS Velocity u rms (cm/s) 


vW - w 


2.4 



with neutrally buoyant glass beads of 10 /im mean diameter. The local velocity of the fluid is 
captured by a high speed camera that records the movement of particles over a square area 
of side length L = 7.7 cm. The laser sheet illuminating the particles is roughly 1 mm thick, 
but the camera's depth of field is a tenth of this. Thus only the horizontal components of 
velocity are captured. The buffer size of the camera limits the duration of data collection to 
5 s at a frame rate of 400 frames per second. The velocity field is constructed by correlating 
pairs of images with a particle tracking program. The spatial resolution of the measurements 
depends upon the length of a pixel in the camera's sensor. Hence the spatial resolution in 
this experiment is 7.7 cm/1024 pixels = 75 /im. The temporal resolution is the sampling rate 
which is 7.5 ms for this experiment. Hence the velocity resolution obtained in this experiment 
is 1 cm/s. To ensure reliability of the measured velocities, every n th and {n + 3) rd image is 
employed in velocity field construction. The weak non-zero mean velocity observed in the 
flow is systematically subtracted out. Physical parameters that characterize the turbulence 
are given in table UJ 

III. RESULTS AND DISCUSSION 

Though the Reynolds number in the experiment is low, the flow exhibits a well-defined 
inertial range. Figure El shows the space and time-averaged S$(r) = (A(wi(r)) 3 ) vs. r in 
the range 0.07 cm < r < 3.5 cm. Each instantaneous third-moment is constructed from an 
average over all particle pairs in the field of view. There are roughly 10 5 particle pairs at 
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each r- value contributing to the spatial average at each instant. Twelve time-uncorrelated 
measurements of the third-moment are averaged to obtain the plot in Fig. El The interval 
between these measurements is 1.4 s, which is the lifetime of the largest eddies, as determined 
from the velocity autocorrelation function. Therefore each data-set provides 4 measurements 
uncorrelated in time. Three such data-sets are taken in quick succession to obtain the 
average. In forming this spatial average of Ss(r), data for r < 0.07 cm are excluded, since 
velocity differences cannot be reliably measured below this value. The statistical error in 
the time averaged S3 (r) is less than 5% as observed from the error bars in fig. 

Shown in Fig. El is x 3 P(x) where x = Aui{r,t) at the three indicated r-values 0.3, 0.5 
cm, and 0.7 cm. The measurements were made at one instant of time. The integral of this 
function over x is the the instantaneous third moment of interest. Each P(Av,l(t, t)) is 
constructed from 10 5 particle pairs. The statistical fluctuations in each of the curves is ade- 
quately small to yield a statistically significant value of the (fluctuating) third moment and 
hence the pdf of this moment P(S , 3(r)), at least for the values of r listed above. According 
to Eq. the pdf of the energy flux is then given by P(e) = P(S , 3(r))/(— ^r). 

The inset to Fig. H] shows a 5-second time record of e(r,t) for r = 1.925 cm. One sees 
that the energy flux is positive most of the time, but does show intervals of reversed energy 
flow. The figure itself shows the pdf of energy flux at the same three inertial range r-values 
of Fig. El One expects that in the inertial range all three curves should coincide, and they 
approximately do so. This r-independence in the inertial range is expected for the third- 
moment. The energy flux here is experimentally obtained from Eq. |2]for a sub-domain B of 
size R = 1.925 cm. As required by Eq. El the pdfs are positively skewed. 

The energy flux is also experimentally obtained for sub-domains B of size R {R= 7.7 cm, 
3.85 cm and 1.925 cm) keeping the particle separation r = 0.7 cm constant. The "width" 
of their pdfs decreases with increasing R, as one would expect. In general agreement with 
previous observations, it is found that the probability of backscattering is significant, but it 
is measured to be a monotonically decreasing function of R when R increases beyond l . 

IV. NUMERICAL SIMULATIONS 

The experimental results have been supplemented by a numerical simulation, using the 
SPH algorithm to simulate flows with particles . In this approach, the fluid is represented 
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by particles, whose positions, and velocities, Vj evolve according to the equation of motion 



§ = Vj ; ^ = _ V p, + ,, + D i (3) 

where Vp^ is the pressure gradient, fj the forcing term, and Dj the dissipation term. The 
pressure is obtained from the local density pi from the equation of state p = p 7 (7 = 9), 
which corresponds to a very weakly compressible gas provided |i>j| is small compared to the 
velocity of sound. The density p as well as the gradients are computed at a point x by 
interpolating the properties of the particles with a kernel function characterized by a size 
h. In this approach, scales smaller than h are unresolved. The SPH method allows one to 
simulate flows with a resolution no better than h. 

The issue is thus to adequately parametrize the sub-grid energy dissipation e h occurin 
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below the scale h. The precise form for Eh is inspired by the Kolmogorov equation 
it corresponds to an energy dissipation of : 

^ = -J^(( v ij- r ij) v D- ( 4 ) 

Here the angular brackets denote an average over all particles and this term is negative 

on the average, and Vjj = (vj — Vj) is the velocity difference between the particles % and 

j separated by the distance r^. The parameter v t is a dimensionless constant deduced to 

be of the order of unity Q], and should not be confused with kinematic viscosity. The 

experimental apparatus permits measurement of the postulated eddy damping term, and a 

comparison with the form of Eq. El Agreement between the flux in Eq. Eland the dissipation 

in Eq. 0]is extremely good. With a properly adjusted value of u t , the correlation between e(t) 

and Eh is ~ 0.98 at the three values of R = 7.7, 3.85 and 1.925 cm. All the properties of the 

energy flux computed with Eq. HI are essentially identical to the properties computed with 

Eq. |21 as shown in Fig. [5]. The precise value of v t is ~ 4.8. The value of v t depends on the 

precise geometry of the subdomain B. This experimental finding is a striking confirmation 

I I 

of the validity of the eddy-damping postulated in [14] . 

The implementation of SPH algorithm is described in Q].In the calculation, the forcing 
term is chosen as a random superposition of three low-wavenumber fourier modes. The 
velocity of each particle was no larger than 1/10 the velocity of sound, so the flow is effectively 
incompressible. The energy flux over a sub-volume of the system is computed numerically, 
as effectively done in the experiment. The total energy flux averaged over a system volume 



of side- length L/2 (for a simulation box of side- length L) has large fluctuations. The inset to 
Fig-Et a ) shows that the fluctuations remain negative for as long as the simulations have been 
carried out. When averaged over a smaller fraction of the volume of side-length L/8, one 
finds that positive fluctuations of the averaged energy-flux become possible. This is evident 
from the inset to Fig. G^b) and the corresponding pdf. The smaller the box size larger 
are the fluctuations. This finding is fully consistent with the experimental observations. 
As in the experiment, it was found that the probability of energy flux backscattering is a 
monotonically decreasing function of the scale R over which the flux is defined. 



V. SUMMARY 



Fluid turbulence research has focussed either on small scale properties, with the objective 
of understanding "intermittency" , or on averages of global properties of interest for most 
applications. Only recently has a more systematic investigation of global fluctuations been 
undertaken The focus here is on the energy transfer to smaller scales, a hallmark 

of hydrodynamic turbulence. Whereas on average, the fluid transfers energy to smaller scales 
via the cascade process, the observed fluctuations are so large, as to actually reverse the 
energy flux from small to large scales This experimental study of the fluctuations of 
the energy flux is based on an experiment carried out at moderate Reynolds number. The 
results are based on Eq. El It has been shown that the probability of energy backscattering 
decreases when the size of the system increases. Also, the form of the eddy-damping term as 
proposed in [14( has been justified. The experimental results are corroborated by the results 
of the calculation based on the particle method. 

The systematic study of scale dependence of the fluid properties is a novel aspect of this 



work. It would be interesting to understand these results in the light o 
theory derived in the context of nonlinear, out-of-equilibrium systems (2 



recent fluctuation 
2( 
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FIG. 1: a) Top view and b) Side view of the experimental setup. A tank of lateral dimensions 
1 m x 1 m is filled with water to a depth of 30 cm. A laser beam passing through a cylindrical 
lens generates a sheet of laser that illuminates a horizontal plane of the turbulent fluid 4 cm below 
the surface. Turbulence is generated by a system of 36 capped rotating jets situated at the tank 
floor. Neutrally buoyant tracer particles are suspended in the turbulent fluid. A high speed camera 
suspended above the tank records the motion of tracers as they scatter light upon entering the 
sheet of illumination. 




FIG. 2: (Color online) The time averaged third-moment of longitudinal velocity differences (Ss(r)) 
as a function of particle separation r (cm). The plot is constructed by averaging over twelve 
time-uncorrelated velocity snapshots. 
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FIG. 3: (Color online) The integrand of the third-moment of longitudinal velocity differences can 
be constructed reliably for an instantaneous snapshot for different spatial separations r (here taken 
for r= 0.3, 0.5 and 0.7 cm). 




FIG. 4: (Color online) The inset shows the time trace of energy flux obtained from Eq. |2J The 
main figure shows the pdf of energy flux obtained for three different spatial separations r = 0.3, 
0.5 and 0.7 cm. 
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FIG. 5: (Color online) (a) shows the time-trace (inset) and the pdf of energy flux for a simulation 
box of side length L/2. The fluctuations are always positive and never change sign, (b) shows 
the time-trace (inset) and the PDf of energy flux for a simulation box of side length L/8. The 
fluctuations switch sign frequently when the box size has diminished. 
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